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Abstract 

A method has been developed for analyzing the 
flow in subsonic axisymmetric inlets at arbitrary 
conditions of freestream velocity, incidence 
angle, and inlet mass flow. An improved version 
of the method is discussed and comparisons of re- 
sults obtained with the original and improved 
methods are given. Comparisons with experiments 
are also presented for several inlet configurations 
and for various conditions of the boundary layer 
from insignificant to separated. The paper dis- 
cusses applications of the method, with several 
examples given for specific cases involving inlets 
for VTOL lift fans and for STOL engine nacelles. 

Introduction 

Many proposed advanced aircraft, whether CTOL, 
STOL or VTOL, require propulsion system inlets to 
operate efficiently over wide ranges of flight 
speed, incidence angle and inlet throat Mach 
numbers (mass flow rates). These requirements 
can be quite severe for a fixed-geometry axisym- 
metrlc inlet. Therefore, considerable research 
and development effort is tequlred for the design 
of such inlets. 

The principal tool in inlet design has been 
wind tunnel experiments with scaled model inlets. 
Wind tunnel testing is both lengthy and expensive. 
To minimize the amount of wind-tunnel testing re- 
quired and to enaure that reasonable geometries 
are tested, a reliable theoretical method of In- 
let analysis is needed. The method should be 
able to calculate the potential and vlacous flow 
in Inlets of arbitrary geometry and combinations 
of operating conditions* 

Such a mat hod lias evolved over the past several 
years at the NASA Lewis Research Center. The 
original motivation for a potential flow analysis 
(rtf. 1) was ths need to design Inlets for an in- 
house VTOL lift fan test program (refa. 2(3). 

Ths method was quite successful at this and waa 
extended to sev v raL other applications (ref. 4). 
When the method w*j applied to STOL inlet designs, 
ths boundary layar recamc significant, and the 
boundary layer or viscous calculations were In- 
corporated Into the system (ref, 5). 

A status report on the resulting method as of 
lata 1973 is given In reference 6. Since that 
paper, the met bed waa Improved and many additional 
applications were made (e.g., refa. 748). Ths 
present paper (essentially an update of ref. 6) 
will cover the elements of the method, a com- 
parison of the improved method with the original 
method, several comparisons with experiment, and 
a discussion of number nt design and analysis 
applications. 
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Nomenclature 

A, B, C combination coefficients, eq (1) 

local skin friction coefficient (ratio 
of wall shear stress to dynamic pres- 
sure at edge of boundary layer) 

D diameter 

M Mach number 

P pressure 

5 surface distance 

U T rotor tip speed 

V velocity 

W Inlet mass flow 

cc inlet incidence angle 

6 boundary layer thickness 

6* boundary layer displacement thickness 

A^ change in rotor incidence angle 

p density 

0 source strength 
Subscripts: 

c control station 

cor corrected for compressibility 

cr critical 

h highlight 

1 Incompressible 

j basic solution 

max maximum value 

ref reference value 

s static conditions 

t total or stagnation conditions 

th throat 

* frsastrsam value 

Superscripts: 

average value 

— > vector quean *ty 
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Method of Solution 
Statement of the Problem 

The basic problem to be solved is to calculate 
the compressible viscous flow in an arbitrary 
axisymnetric inlet at any combination of oper- 
ating conditions of inlet mass flow rate, ft, 
freestream velocity V^, and inlet incidence 
angle, a (fig. 1). At non-zero incidence angle 
the flow in and around the inlet is three- 
dimensional. At the present time there is no 
exact practical compressible viscous flow method 
of solution (computer program) capable of handling 
this inlet problem. Therefore, the problem is 
solved in several steps (shown schematically in 
fig. 2) as follows: 

1. Geometry representation (Program SC1RCL) 

2. Incotroresslble potential flow basic solu- 
tions (EOD) 

3. Combined solutions with compressibility 
correction (COMBYN) 

4. Boundary layer calculations (VISCUS) 

5. Iterative loop 

The four computer programs are available from 
COSMIC , Computer Center, Information Services, 

112 Barrow Hell, University of Georgia, Athens, 
Georgia, 30602. Programs SC1RCL, EOD and COMBYN 
are one unit (ref. 9) with number LEW-12152; pro- 
gram VI SCUS (ref. 5) is number LEV-12178. Each 
■tap will be next described in some detail. 

Ceomtn/ Representation 

The inlet le assumed to be axlsynsetrlc and is 
represented by its meridional profile. This pro- 
file is broken Into segments at convenient tangent 
points as shown in figure 3. Each segment may be 
defined by an analytic expression or a set of 
points. The Inlet duct * wells and the outer sur- 
face (nacelle or bellmouth) must be extended far 
downstream (fig. 3) to facilitate obtaining 
accurate potential flow solutions in the inlet 
region of Interest. The geometry program (SC1RCL, 
fig. 2) prepares coordinate-point Input for 
efficient use of the potential flow program and 
also prints out Information such as curvature, 
well angles, flow area distribution, etc*, which 
Is useful in preliminary screening of proposed 
inlet shapes. 

Incompressible Potential Flow Basic Solutions 

The Douglas-Net imann program (rsf. 10) Is used 
for calculating the Incompressible potential flow. 
Briefly, the program util lass e distribution of 
sources or sinks of initially unknown strength o 
to represent the inlet profile. This represen- 
tation results in an lntsgrsl aquation (sea rsf. 

10 for datalls) which Is sxsct for s continuous 
distribution of sourcs strength. This continuous 
distribution la approximated by representing the 
inlet profile by a finite number of discrete 
elements characterised by s point on the element 
(e»|*, the midpoint) called the control point, 
lech element has the asms predetermined type of 


source strength distribution (e.g., constant, 
linear, parabolic). This approximation results in 
a set of linear algebraic equations that are solved 
by matrix metnods for the source strength at the 
control points. Velocities at the control points 
and at specified off-body points are then calcu- 
lated from the source distribution. 

Method of approximation . Two methods of approxi- 
mation have been used as Bhovn in figure 4: (1) 

the original method, called the base method, which 
has been in use at NASA Lewis for several years; 
and (2) the improved method, called the higher 
order method which was recently put into use. 

The base method (ref. 10) uses flat (linear) 
surface elements and assumes constant, source 
strength over each element (fig. 4(a)). To obtain 
solutions of adequate accuracy this method often 
requires very large numbers of elements and con- 
sequent long computer times. 

The higher order method (ref. 11) uses curved 
(parabolic) surface elements and assumes a linear 
variation in source strength over each element 
(fig. 4(b)). For a given accuracy this method 
requires fewer elements than the base method with 
consequent savings in computer time. Conversely, 
a greater accuracy can be obtained within the 
element -number limitations of a given program- 
computer system with this method. 

Types of basic solutions . The Douglas potential 
flow program is used to obtain three basic solu- 
tions which are used in linear combination (to 
be explained under "Combined Solution" below) In 
order to satisfy the prescribed operating condi- 
tions (fig. 1). Two types of sets of basic 
solutions, as shown In figure % have been used at 
NASA Lewis. The first is the closed-duct meLhod 
(In uso for several years) and the second is the 
shroud-vortlcity method (recently put into use). 

The closed duct method uses a combination of a 
closed-duct inlet (fig. 5(a)) and an open-duct 
inlet both in an axial freestream flow to obtain 
a static arbitrary maas flow. Thla method has 
some shortcomings that will be illustrated later 
under "Comparison of Methods of Solution". 

The ahroud-vorticity method fflg. 5(b)) utilises 
a distribution of unit vorticas (in addition to 
tha distribution of sources that repreaant the 
inlet profile) on the shroud surface to i:*duce a 
static mass flow through the inlet. Any arbitrary 
static mass flow can be obtalnad by the use of a 
multiplicative factor. This method does not suffer 
tha shortcomings of the closed-duct method. 

Possible solution procedures . Tha two methods 
of approximation and tha two me thole of inducing 
static maas flow are independent ot eech other so 
that whan they are used In all possible combin- 
ations there results four different procedures 
for tha solution: 

1. Higher order approximation * shroud 
vortlclty. 

2. Higher order - closed duct 

3. Baas method - shroud vortlclty 
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4, Base method - closed duct 

The effects of these four procedures are inves- 
tigated in reference 12. In that reference the 
four methods are applied to both STOL and VTOL 
inlets at several operating conditions. The 
greatest difference In results (static pressure 
distribution) are obtained between method 1 (the 
method currently In use) and method 4 (original 
method) for thin-walled nacelle type Inlets in 
static operation. An example is shown in figure 
6. The differences are quite significant, with 
a particularly anomolous behavior for the old 
method in the region of the fan face. The loca- 
tion of the anomoly Is the belnnnlng of the 
minimum thickness of the nacelle* The element 
slse here Is too large relative to the nacelle 
thickness* The anomolous behavior could be elimin- 
ated, and the curve of method 4 in figure 6 
could be made to agree with that of method 1 by 
using a sufficiently large number of elements. 
However, for this particular inlet that number 
exceeds the limitations of the program. The con- 
clusion from the Investigation of reference 12, of 
which figure 6 is an example, la that the new 
method gives greeter accuracy with fewer elements 
and consequent smaller computing time* 

Combined Solution 

The incompressible potential flow computer pro- 
gram Just discussed is used to obtain three in- 
dependent basic solutions (fig. 5). These three 
basic solutions V., j - 1,2,3 are then combined 
(program COMBYM, fig. 2) into a solution of Inter- 
est v having arbitrary flow conditions of V w , a, 
and mess flow W (fig. 1). The combination equation 
is: 

^ ( 1 ) 

where A, B, end C are obtained t.*om the three flow 
conditions. Thus, once the basic i'ow solutions 
are obtained for a specified geometiy, any solu- 
tion of interest for that geometry jan be obtained 
without repeating the more time-consuming potential 
flow calculations* 

Compressibility 


gram VISCUS, fig. 2). Reference 5 contains a coi&- 
plete documentation of program VISCUS and refer- 
ences to the original sources. Program VISCUS 
calculates jjoundary layer profiles, displacement 
thickness 6 , skin friction coefficient Cj , etc., 
at each station, and also predicts transition 
from laminar to turbulent flow. Separation 
(whether laminar or turbulent), is predicted when 
Cf is zero. 

The boundary layer calculations are based on 
the assumption that the flow la two-dimensional, 
either planar or axlsymmetrlc. Questions have 
been raised (e.g., by Schaub and by Presley and 
Perkins, DISCUSSION, of ref. 6) about the validity 
of using two-dimensional boundary layer calcu- 
lations for flow that is three-dimensional. Up 
till now, the boundary layer calculations were 
made only along the inlet windward meridian (lower 
lip in fig. 1) because the longitudinal velocity 
gradients are most severe there and separation 
will occur there first. An attempt to assess the 
validity of using the two-dlmenslonsl boundary 
layer calculation along the windward and other 
meridians by comparing the calculated boundary 
profiles with experimental at several circumfer- 
ential locations is currently nearing completion 
(ref. 14). In another approach, a study is under- 
way to calculate the boundary layer along three- 
dimensional streamlines calculated from the 
potential flow on the inlet surface. This approach 
should result in the best approxlmstlon possible 
with the existing program system. Profiles calcu- 
lst ed by this method will be compared with experi- 
mental data to assess the usefulness of this 
approach. 

Another shortcoming of the boundary layer cal- 
culation is the neglect of shock Interactions* 

Since many cases of current Interest contain 
rsglona of local supersonic flow, it may be neces- 
sary to account for possible shock-boundary- layer 
interaction In the boundary layer calculations. 
Furthermore, the transition model is not able to 
predict separation bubbles that appear to be pre- 
sent in the experimental data sspsc tally in small- 
scale. However, It might be feaelble to Include 
the predictions of s separation bubble in as im- 
proved transition model. 


The velocity obtained by equation (1) la incon- 
presslble and is corrected for compressibility by 
the Liable in-Stockman compressibility correction 
(ref* 13). 


cor 


(fc) 


v 


( 2 ) 


where ell the terme on the right hand elde are 
obtained from the lncompreeelble flow eolation 
or the Input flew conditions. This correction 
require# no alteration of the inlet geo net '7 -rod 
It can handle local sonic and supersonic velo- 
cities. From the compressible velocity, V , 
other flow properties (Mach number, pressure 
retie, streamlines, etc.) ere obtained. 


The surface Mach number distributions obtained 
f ran program CONST* are used ee Input to the 
Herrlng-Meller boundary layer calculation (pio- 


Itcratlve Loop 

If the boundary layer le significant in the in- 
lst rsglon of intsrsst. It nay be deelreable to 
add the dleplecenent thickness 6 * to the original 
Inlet prof 11 s end repeat the entire solution pro- 
cedure, thus obtaining e new Mech numbs r distri- 
bution, new £*, etc* This process nay be iterated 
to satisfactory convergence. Usually one iteration 
is sufficient. In parametric studies or prelimi- 
nary design screening, often no Iteration le needed 


To indicate the accuracy of the various aspects 
of the prediction net hod several comparisons with 
experimental date will be given. These will range 
over Incompressible flow, compressible flow with 
insignificant boundary layer, end compressible 
viscous flow with various conditions of the 
boundary layer. 
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Incompressible 

The surface pressure distributions In » chord- 
wise cut of s VTOL fan-ln-wing inlet st low 
forward velocity are shown In figure 7. The 
experlaiental data of figure 7 were obtained from 
the National Research Council of Canada (ref. 15). 
Three surfaces are shown on the plot: the forward 

surface of the bsllmouth; the centerbody; and the 
aft surface of the bellmouth. The agreement 
between theory and experiment la quite good every- 
where on the Inlet. This case la Included to 
Illustrate the adequacy of the method of geometry 
representation even when the calculatlonal model 
(see fig. 3(b)) dlffars rather significantly from 
tha real geometry. 

Compressible 

The next ceee illustrates the applicability of 
the method when the flow is compressible. Figure 
8 shows the theoretical and experimental surface 
pressure distributions oi. a VTOL fen-in-pod In- 
let (ref. 3). Both Incompressible snd incompres- 
elble-correctad-for-compressibllity theoretical 
curves are given. The experimental static pres- 
sures agree quite well with the theory corrected 
for compressibility along the entire surface of 
the Inlet. 

This good agreement when the flow is compres- 
sible Is usually obtained with other configurations 
both VTOL and STOL, and at various operating con- 
ditions provided thsrs are no extended regions of 
local supersonic flow. That there may be dis- 
agreement when there le supersonic flow it not 
surprising since the relation between local 
etreamtube area and velocity is different for 
supersonic ans subsonic flows and the method Is 
based on the subsonic relation. Furthermore, 
there may be other real flow effects present In 
a supersonic flow not accounted for in the ana- 
lysis such as the presence of shocks or shock- 
boundary- layer int erect ion. 

To illustrate the kind of agreement that can 
be expected when regions of local supersonic flow 
are present, two examples will be given of STOL 
engine inlets. First e model STOL engine inlet 
st approach operating condition la shown in 
figure 9. The figure shows a comparison of cal- 
culated end measured static prtssurea on both 
the windward and lssward sides of the inlet. 

The agreement would probably be quite adequate 
for most design or analysis applications. 

Another small scale model STOL engine inlet 
operating at a takeoff throat Mach number of 0.76 
(for noise reduction purposes) in t 34-knot 
crosswind (o - 90°) is shown in figure 10. Here, 
although the theory follows the data in a general 
way, there are two regions of disagreement: (1) on 
the lip; and (2) in the rearward portion of the 
diffuser. The disagr eem e n t in the diffuser le 
due tc neglect lag the boundary layer thickness in 
the calculations (as will be seen in the next 
section). The disa gr eeme n t in the region of the 
internal lip may bo due to the presence of a 
separation bubble near the highlight and a shock 
(or compression waves) near the throat . 


It should also be remembered that the observed 
disagreements between theory and small-scale wind 
tunnel experiments are not necessarily the fault 
of the theory. There may be errors In the mea- 
surements. model contours, wind tunnel wall effects 
or other induced effects. For example, in some 
cases of poor agreement between theory and experi- 
ment, the agreement could be improved markedly by 
inputlng into the calculation an incidence angle 
two or three degrees different from the measured 
value, thus indicating a possible induced effect 
in the wind tunnel tests. In fact, a alight in- 
crease in incidence angle would improve the agree- 
ment on both windward and leeward sides of the 
inlet of figure 9. 

Viscous (Boundary Layer) 

A small-scale translating-centerbody sonic Inlet 
in the retracted (cruise) configuration (fig. 11(a)) 
is used for comparing theoretical results including 
boundary -layer effects to experimental data. Var- 
ious conditions of the boundary layer are obtained 
by varying the inlet incidence angle a of the 
model. In all cases, the one-dimens ionel throat 
Mach number is 0. 50 and the freestream Mach number 
is 0.13 (These results are all taken from ref. 8.) 

Boundary layer attached . Surface Mach number 
distributions for zero incidence angle are shown 
in figure 11(b). Two theoretical curves are 
shown to Illustrate the effect of the calculated 
boundary layer on the surface Mach number. The 
solid curve is the potential flow (no 6* correc- 
tion); the broken curve Is "he potential flow 
obtained with < 5 * added to t e inlet profile. With 
6 * added, there la an excellent agreement with the 
theory in the diffuser region. No separation is 
Indicated either by theory or by experiment. 

Diffuser separation . Results are shown in 
figure 11(c) for an incidence angle of 40°. Here 
the theory, with or without the correction, 
predicts separation in the diffuser at about the 
came point as the experimental data indicate the 
start of a "pressur -lateau” commonly associated 
with an actual sep*.t«tion. 

Lip separation . On figure 11(d) results are 
given for an inlet incidence angle of 50°. The 
theoretical Mach number distribution is what would 
be obtained If the inlet flow did not separete. 

The experimental data clearly indicate eeperation 
on the Inlet lip. The theory also predicts lip 
separation, and although poeelbly not at the exact 
point ae the actual separation, it la quite ade- 
quate ae a guide In Inlet design. 

In canes where the theory predicts separation 
the calculations atop and thare la no 6* calcu- 
lated beyond the point of separation. For caaat 
whore a 6* correction la used, the distribution 
of £* Into the separated region can be obtained by 
extrapolation using an unseparated case (lower 
a) as a gulda (see raf. 8). 

Additional comparison for this and othtr con- 
figurations can ba found in rsfarence 6. On tha 
basis of these results, the theory seems adequate 
in predicting the boundary layer behavior for this 
type of Inlet configuration. 
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Applications 

The calculation method has been used In many 
applications. Among them are: design of inlets 

for test models; analysis of existing inlets; 
parametric studies to aid in future design work; 
and test support. Examples of test support are 
calibration of inlets for mass flow rate measure- 
ment and interpretation of experimental results 
by comparison with the calculated results. 

In this paper several applications of the cal- 
culation method will now be discussed, some 
utilizing only the potential flow part of the 
calculation and others using also the boundary 
layer calculations. 


the potential flow due to the presence of the 
rotor.) The effect of potential flow inflow dis- 
tortion on fan stage performance is illustrated In 
reference 2. Additional sets of calculated inci- 
dence contours can be obtained to study the effect 
of different design parameters such as inlet depth, 
transition velocity, inlet profile, and rotor tip 
speed. 

STOL Inlets . The next two applications are taken 
from a theoretical screening study of Inlets for 
the Quiet Clean Short-Haul Experimental Engine 
(QCSEE) Project (ref. 7). The first is a study of 
the effect of varying internal-lip contraction 
ratios and the second is a study of the effect of 
external forebody shape on the Internal flow. 


Potential Flow 


VTOL inlets . Two applications of the method to 
VTOL inlets will be discussed: one a design 

application; the other an analysis application. 

An example of the use of surface velocity dis- 
tributions and passage velocity profiles in the 
design of a VTOL fan-ln-pod inlet la shown in 
fig. 12. The specific problem in this case is 
the determination of a near-optimum location of 
the point of tangency between the bellmouth and 
the pod surface. Three locations of this tangent 
point are shown In the inset in fig. 12(a). Tht 
theoretical surface velocity distributions on 
the three different bellaouths are also shown for 
both static and crossflow operation. It can be 
seen that both the velocity peaks and the unfavor- 
able velocity gradients in crossflow are reduced 
as the tsngsnt point Is moved out to a larger 
radius. Thus case C would be expected to have the 
best crossflow performance. However, at static 
conditions case C shows s hlghsr velocity peek 
end s mors adverse velocity gradient than cases 
A and B. In ordsr to select a beat shape a com- 
promise may be made between the static and cross- 
flow operation. 

The radial velocity profiles at the fan face 
are ahovn in fig. 12(b) for both static and cross- 
flow conditions. The differences between the three 
cesee were not significant enough to affect the 
choice. However, another application of the 
method can be pointed out here, namely that the 
calculated static velocity profile as well ee the 
upstream streamlines can be used as input to the 
fan rotor design. 


A second example (fig. 13) la the calculation 
of VTOL rotor Inflow distortion In crossflow for 
the fan- in-wing configuration shown. The rotor 
of the fan and the Inlet were both deelgned for 
otatlc operation with a ratio oi tip speed U T to 
fen axial velocity V c of 1.7. U this Inlet is 
operated at a ratio of crossflow velocity to fan 
axial velocity of around 0*4, the Incidence angle 
of the potential flow relative to the rotor blades 
will deviate from the design value by magnitude A. , 
as indicated by tbs contours in figure 13* It can 
be soon that. In the plana of the rotor Inlet, the 
incidence eagle dietortlon due to the potential 
flow alone can be severe, (The incidence angle 
distortion ee calculated does not Include inlet 
total praaaure variations or the modification of 


The first example study was undertaken to ensure 
that a reasonable range of contraction ratios was 
chosen for experimental evaluation. Figure 14 
shows the effect of the internal lip contraction 
ratio, (Dy | /D t - )^, on the surface Mach number dis- 
tribution at qCSEE takeoff conditions (Mj h * 0.79, 
M*, ■ 0.12, a - 50°). The figure shows that the 
higher tne contraction ratio, the better the Mach 
number distribution, l.e., the lower the peek value 
and the less severe the unfavorable gradient. 
(However, the takeoff performance must be compro- 
mised with the cruise performance, and the highest 
contraction ratio lip may not be best overall.) 


The second example taken from the QCSEE s*ud> 
of reference 7 is the affect of external forebody 
bluntness on the peak Mach number in the internal 
lip surface (fig. 15). Two of the internal lips 
of figure 14 (with contraction ratios of 1.46 end 
1.56) were studied with three different external 
forebodies. The external forsbodles are charac- 
terized herein only by their value of bluntness. 
Bluntness is a parameter which takes into account 
the local bluntness of the forebody near the high- 
light, the aspect ratio of tha forebody, and the 
relative scale of the forebody thickness to the 
highlight radius. (For details of the forebody 
contour and the quantitative definition of bluntneos 
aee ref. 7) These combinations were analyzed at 
two different operating conditions as shown on 
figure 14. At the QCSEE takeoff conditions (fig. 
15(a)), the peak Mach numbar is quits sensitive 
to external forebody bluntneas with the lower con- 
traction ratio (1.46) lip being more sensitive 
then the 1.56 lip. However, et conditions more 
representative of CTOL takeoff, l.e., an incidence 
angle of 30° and a freestream Mach number of 0.16 
(fig. 15(b)), the internal lip peek Maen number is 
essentially Independent of external forebody blunt- 
ness. 


Viscous Flow 


The two principal uees of the boundary layer cal- 
culations in inlet analysis have been to obtain a 
more realistic potential flow solution by accounting 
for the displacement thickness, and to predict 
boundary layer separation. These usee have already 
been Illustrated la the comparison with experiment 
In figure 11* The criterion for boundary layer 
separation In the method la that tbs local akin 
friction coefficient C* go to aero. Two uees of 
tha akin friction coefficient distribution will be 
discussed In thle section! one Illustrating tha 
affect of Incidence eagle on distribution of Cf ; 
tha other Illustrating tha affect of seals • 
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Effect of Incidence angle . Figure 16 (taken 
from ref, 8) shows the effect of varying incidence 
angle on the skin friction distribution of the 
inlet of figure 11. A local minimum with a low 
value of Cf is a point of likely separation as 
the operating conditions become more severe. The 
zero Incidence curve Indicates two regions of 
possible boundary layer separation (at the two 
minimum*), one in the diffuser (where the flow is 
turbulent) and the other on the internal lip (where 
the flow is initially laminar). (The rise In C f 
from the first minimum Is due to the transition 
from laminar to turbulent flow.) 

At an Incidence angle of 20°, separation Is In- 
dicated in the diffuser. As Incidence angle Is 
increased further up to 40° the separation point 
moves upstream somewhat but remains in the diffuser 
at a position close to the diffuser wall inflec- 
tion point (S * 0*6). At the same time, the mini- 
mum on the lip decreases toward zero. Eventually, 
at 50° incidence, the separation point occurs at 
the lip as the minimum reaches zero there. That 
lip separation Is predicted by the first minimum 
going to zero rather than the second minimum moving 
forward has been verified by taking very small in- 
crements in Incidence angle (ref. 16). In refer- 
ence 16 It was found that the forward movement of 
the diffuser separation location eventually stops 
and furthar increases In incidence angle reduce 
the lip minimum until It reaches zero Indicating 
separation. 

Experimental results correlstlng total pressure 
loss with sspsrstlon location (ref. 8) for this 
Inlet indicate that diffuser separation in the 
region of S • 0.6 leads to only small losses, 
whereas lip separation leads to Intolerable losses. 
Thus, It appears that a distribution of Cf, like 
that of figure 16, Is undeslreable and If possible 
the inlet design should be modified to raise or 
eliminate the minimum on the lip. 

Scale affects . The skin friction distribution 
is used in reference 17 to investigate the effects 
of scale on boundary layer separation prediction. 
Figure 17 (taken from ref. 17) shows ths skin 
friction distribution for ths inlet geometry of 
figure 10 at several different scales. The smaller 
ths seals, ths mors likely is the boundary layer 
to separate. If ths seals were mads sufficiently 
small, ths separation point would occur at the 
lip as In the esse cf increasing angle of attack 
(fig. 16). 

Ths seals affects shown in figure 1/ do not In- 
clude the effects of seals on possible shock- 
boundary layer Interaction or on more complex 
t renal t ion mechanisms (such aa laminar separation 
with turbulent reettachment), sines these affects 
are not currently included in ths boundary layer 
program. The effects of longitudinal curvature 
on transition or boundary layer eddy viscosity 
era also not included in these particular calcu- 
lations. 

Figure 17 and othnr result a of reference 17 in- 
dicate that tasting at small scale tends to give 
pessimistic estimates of boundary layer behavior. 

At attempt will be mads in reference 17 to provide 
a proc e d u r e for predicting b ou n d ary layer at full 


scale based on a combination of small scale test 
results and an improved version of the boundary 
layer calculations. 

Concluding Remarks 

A method for the prediction of potential and 
viscous flow in subsc ilc axlsymmetrlc propulsion 
system Inlets at arbitrary Incidence angle and flow 
conditions has been developed. The method has 
already proven to be a very useful and powerful 
tool for both analysis and design purposes. The 
various elements of the method are constantly being 
updated, thus improving both its versatility and 
accuracy. It is probably the best tool available 
at present for analyzing the compressible viscous 
flow In axlsymmetrlc subsonic Inlets at arbitrary 
incidence angle and prob* K ly will continue to be 
until exact three-dimens K.nal compressible Inlet 
flow computer programs become practical. 
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Figure L - Inlet geometry and flow conditions: inlet mess fleer rete, 
W; free stream velocity, and inlet Incidence angle, o . 
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Figure 6. - Effect of solution procedure on pressure distr button 
on a STOL engine inlet Static condition M, h -0.77. 
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Figure 8. - Comparison with experiment of fan -In -pod inlet 
surface pressure distribution. Static condition, V* - a’ 
Average rotor inlet Mach number, A c ■ 0.57. 







Figure 1L - Inlet with boundary tayer, 0.5(1 
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